Search for: All records

Abstract The metallicity of galaxies, and its variation with galactocentric radius, provides key insights into the formation histories of galaxies and the physical processes driving their evolution. In this work, we analyze the radial metallicity gradients of star-forming galaxies in the EAGLE, Illustris, IllustrisTNG, and SIMBA cosmological simulations across broad mass (10^8.0M_⊙≤M_⋆ ≲ 10^12.0M_⊙) and redshift (0 ≤z≤ 8) ranges. We find that all simulations predict strong negative (i.e., radially decreasing) metallicity gradients at early cosmic times, likely due to their similar treatments of relatively smooth stellar feedback not providing sufficient mixing to quickly flatten gradients. The strongest redshift evolution occurs in galaxies with stellar masses of 10^10.0–10^11.0M_⊙, while galaxies with stellar mass < 10¹⁰M_⊙and >10¹¹M_⊙exhibit weaker redshift evolution. Our result of negative gradients at high redshift contrast with the many positive and flat gradients in the 1 < z < 4 observational literature. Atz > 6, the negative gradients observed with JWST and the Atacama Large Millimeter/submillimeter Array are flatter than those in simulations, albeit with closer agreement than at lower redshift. Overall, we suggest that these smooth stellar feedback galaxy simulations may not sufficiently mix their metal content radially, and that either stronger stellar feedback or additional subgrid turbulent metal diffusion models may be required to better reproduce observed metallicity gradients. 

Abstract We present 0.6–3.2 pc resolution mid-infrared (MIR) JWST images at 7.7μm (F770W) and 21μm (F2100W) covering the main star-forming regions of two of the closest star-forming low-metallicity dwarf galaxies, NGC 6822 and Wolf–Lundmark–Melotte (WLM). The images of NGC 6822 reveal filaments, edge-brightened bubbles, diffuse emission, and a plethora of point sources. By contrast, most of the MIR emission in WLM is pointlike, with a small amount of extended emission. Compared to solar-metallicity galaxies, the ratio of 7.7μm intensity ( $I_{\nu }^{\mathrm{F770W}}$), tracing polycyclic aromatic hydrocarbons (PAHs), to 21μm intensity ( $I_{\nu }^{\mathrm{F2100W}}$), tracing small, warm dust grain emission, is suppressed in these low-metallicity dwarfs. Using Atacama Large Millimeter/submillimeter Array CO(2–1) observations, we find that detected CO intensity versus $I_{\nu }^{\mathrm{F770W}}$at ≈2 pc resolution in dwarfs follows a similar relationship to that at solar metallicity and lower resolution, while the CO versus $I_{\nu }^{\mathrm{F2100W}}$relationship in dwarfs lies significantly below that derived from solar-metallicity galaxies at lower resolution, suggesting more pronounced destruction of CO molecules at low metallicity. Finally, adding in Local Group L-Band Survey 21 cm Hiobservations from the Very Large Array, we find that $I_{\nu }^{\mathrm{F2100W}}$and $I_{\nu }^{\mathrm{F770W}}$versus total gas ratios are suppressed in NGC 6822 and WLM compared to solar-metallicity galaxies. In agreement with dust models, the level of suppression appears to be at least partly accounted for by the reduced galaxy-averaged dust-to-gas and PAH-to-dust mass ratios in the dwarfs. Remaining differences are likely due to spatial variations in dust model parameters, which should be an exciting direction for future work in local dwarf galaxies. 

The property of star formation rate (SFR) is tightly connected to the amount of dense gas in molecular clouds. However, it is not fully understood how the relationship between dense molecular gas and star formation varies within galaxies and in different morphological environments. Most previous studies have typically been limited to kiloparsec-scale resolution such that different environments could not be resolved. In this work, we present new ALMA observations of HCN(1−0) at 260 pc scale to test how the amount of dense gas and its ability to form stars varies with environmental properties. Combined with existing CO(2−1) observations from ALMA and Hαfrom MUSE, we measured the HCN/CO line ratio, a proxy for the dense gas fraction, and SFR/HCN, a proxy for the star formation efficiency of the dense gas. We find a systematic > 1 dex increase (decreases) of HCN/CO (SFR/HCN) towards the centre of the galaxy, and roughly flat trends of these ratios (average variations < 0.3 dex) throughout the disc. While spiral arms, interarm regions, and bar ends show similar HCN/CO and SFR/HCN, on the bar, there is a significantly lower SFR/HCN at a similar HCN/CO. The strong environmental influence on dense gas and star formation in the centre of NGC 4321, suggests either that clouds couple strongly to the surrounding pressure or that HCN emission traces more of the bulk molecular gas that is less efficiently converted into stars. Across the disc, where the ISM pressure is typically low, SFR/HCN is more constant, indicating a decoupling of the clouds from their surrounding environment. The low SFR/HCN on the bar suggests that gas dynamics (e.g. shear and streaming motions) can have a large effect on the efficiency with which dense gas is converted into stars. In addition, we show that HCN/CO is a good predictor of the mean molecular gas surface density at 260 pc scales across environments and physical conditions. 

ABSTRACT The scatter about the mass-metallicity relation (MZR) has a correlation with the star formation rate (SFR) of galaxies. The lack of evidence of evolution in correlated scatter at z ≲ 2.5 leads many to refer to the relationship between mass, metallicity, and SFR as the Fundamental Metallicity Relation (FMR). Yet, recent high-redshift (z > 3) JWST observations have challenged the fundamental (i.e. redshift-invariant) nature of the FMR. In this work, we show that the cosmological simulations Illustris, IllustrisTNG, and Evolution and Assembly of GaLaxies and their Environment (EAGLE) all predict MZRs that exhibit scatter with a secondary dependence on SFR up to z = 8. We introduce the concept of a ‘strong’ FMR, where the strength of correlated scatter does not evolve with time, and a ‘weak’ FMR, where there is some time evolution. We find that each simulation analysed has a statistically significant weak FMR – there is non-negligible evolution in the strength of the correlation with SFR. Furthermore, we show that the scatter is reduced an additional ∼10–40 per cent at z ≳ 3 when using a weak FMR, compared to assuming a strong FMR. These results highlight the importance of avoiding coarse redshift binning when assessing the FMR. 

ABSTRACT The metal content of galaxies provides a window into their formation in the full context of the cosmic baryon cycle. In this study, we examine the relationship between stellar mass and stellar metallicity (MZ*R) in the hydrodynamic simulations Illustris, TNG, and EAGLE (Evolution and Assembly of GaLaxies and their Environment) to understand the global properties of stellar metallicities within the feedback paradigm employed by these simulations. Interestingly, we observe significant variations in the overall normalization and redshift evolution of the MZ*R across the three simulations. However, all simulations consistently demonstrate a tertiary dependence on the specific star formation rate (sSFR) of galaxies. This finding parallels the relationship seen in both simulations and observations between stellar mass, gas-phase metallicity, and some proxy of galaxy gas content (e.g. SFR, gas fraction, and atomic gas mass). Since we find this correlation exists in all three simulations, each employing a subgrid treatment of the dense, star-forming interstellar medium (ISM) to simulate smooth stellar feedback, we interpret this result as a fairly general feature of simulations of this kind. Furthermore, with a toy analytic model, we propose that the tertiary correlation in the stellar component is sensitive to the extent of the ‘burstiness’ of feedback within galaxies. 

Abstract We present the largest catalog to date of star clusters and compact associations in nearby galaxies. We have performed aV-band-selected census of clusters across the 38 spiral galaxies of the PHANGS–Hubble Space Telescope (HST) Treasury Survey, and measured integrated, aperture-corrected near-ultraviolet-U-B-V-Iphotometry. This work has resulted in uniform catalogs that contain ∼20,000 clusters and compact associations, which have passed human inspection and morphological classification, and a larger sample of ∼100,000 classified by neural network models. Here, we report on the observed properties of these samples, and demonstrate that tremendous insight can be gained from just the observed properties of clusters, even in the absence of their transformation into physical quantities. In particular, we show the utility of the UBVI color–color diagram, and the three principal features revealed by the PHANGS-HST cluster sample: the young cluster locus, the middle-age plume, and the old globular cluster clump. We present an atlas of maps of the 2D spatial distribution of clusters and compact associations in the context of the molecular clouds from PHANGS–Atacama Large Millimeter/submillimeter Array. We explore new ways of understanding this large data set in a multiscale context by bringing together once-separate techniques for the characterization of clusters (color–color diagrams and spatial distributions) and their parent galaxies (galaxy morphology and location relative to the galaxy main sequence). A companion paper presents the physical properties: ages, masses, and dust reddenings derived using improved spectral energy distribution fitting techniques. 

ABSTRACT We present radial gas-phase metallicity profiles, gradients, and break radii at redshift z = 0–3 from the TNG50-1 star-forming galaxy population. These metallicity profiles are characterized by an emphasis on identifying the steep inner gradient and flat outer gradient. From this, the break radius, Rbreak, is defined as the region where the transition occurs. We observe the break radius having a positive trend with mass that weakens with redshift. When normalized by the stellar half-mass radius, the break radius has a weaker relation with both mass and redshift. To test if our results are dependent on the resolution or adopted physics of TNG50-1, the same analysis is performed in TNG50-2 and Illustris-1. We find general agreement between each of the simulations in their qualitative trends; however, the adopted physics between TNG and Illustris differ and therefore the breaks, normalized by galaxy size, deviate by a factor of ∼2. In order to understand where the break comes from, we define two relevant time-scales: an enrichment time-scale and a radial gas mixing time-scale. We find that Rbreak occurs where the gas mixing time-scale is ∼10 times as long as the enrichment time-scale in all three simulation runs, with some weak mass and redshift dependence. This implies that galactic discs can be thought of in two-parts: a star-forming inner disc with a steep gradient and a mixing-dominated outer disc with a flat gradient, with the break radius marking the region of transition between them. 

Abstract We use 0.1″ observations from the Atacama Large Millimeter Array (ALMA), Hubble Space Telescope (HST), and JWST to study young massive clusters (YMCs) in their embedded “infant” phase across the central starburst ring in NGC 3351. Our new ALMA data reveal 18 bright and compact (sub-)millimeter continuum sources, of which 8 have counterparts in JWST images and only 6 have counterparts in HST images. Based on the ALMA continuum and molecular line data, as well as ancillary measurements for the HST and JWST counterparts, we identify 14 sources as infant star clusters with high stellar and/or gas masses (∼10⁵M_⊙), small radii (≲ 5 pc), large escape velocities (6–10 km s⁻¹), and short freefall times (0.5–1 Myr). Their multiwavelength properties motivate us to divide them into four categories, likely corresponding to four evolutionary stages from starless clumps to exposed Hiiregion–cluster complexes. Leveraging age estimates for HST-identified clusters in the same region, we infer an evolutionary timeline, ranging from ∼1–2 Myr before cluster formation as starless clumps, to ∼4–6 Myr after as exposed Hiiregion–cluster complexes. Finally, we show that the YMCs make up a substantial fraction of recent star formation across the ring, exhibit a nonuniform azimuthal distribution without a very coherent evolutionary trend along the ring, and are capable of driving large-scale gas outflows. 

Abstract The CO-to-H 2 conversion factor ( α CO ) is central to measuring the amount and properties of molecular gas. It is known to vary with environmental conditions, and previous studies have revealed lower α CO in the centers of some barred galaxies on kiloparsec scales. To unveil the physical drivers of such variations, we obtained Atacama Large Millimeter/submillimeter Array bands (3), (6), and (7) observations toward the inner ∼2 kpc of NGC 3627 and NGC 4321 tracing 12 CO, 13 CO, and C 18 O lines on ∼100 pc scales. Our multiline modeling and Bayesian likelihood analysis of these data sets reveal variations of molecular gas density, temperature, optical depth, and velocity dispersion, which are among the key drivers of α CO . The central 300 pc nuclei in both galaxies show strong enhancement of temperature T k ≳ 100 K and density n H 2 > 10 3 cm −3 . Assuming a CO-to-H 2 abundance of 3 × 10 −4 , we derive 4–15 times lower α CO than the Galactic value across our maps, which agrees well with previous kiloparsec-scale measurements. Combining the results with our previous work on NGC 3351, we find a strong correlation of α CO with low- J 12 CO optical depths ( τ CO ), as well as an anticorrelation with T k . The τ CO correlation explains most of the α CO variation in the three galaxy centers, whereas changes in T k influence α CO to second order. Overall, the observed line width and 12 CO/ 13 CO 2–1 line ratio correlate with τ CO variation in these centers, and thus they are useful observational indicators for α CO variation. We also test current simulation-based α CO prescriptions and find a systematic overprediction, which likely originates from the mismatch of gas conditions between our data and the simulations. 

Abstract Determining how the galactic environment, especially the high gas densities and complex dynamics in bar-fed galaxy centers, alters the star formation efficiency (SFE) of molecular gas is critical to understanding galaxy evolution. However, these same physical or dynamical effects also alter the emissivity properties of CO, leading to variations in the CO-to-H₂conversion factor (α_CO) that impact the assessment of the gas column densities and thus of the SFE. To address such issues, we investigate the dependence ofα_COon the local CO velocity dispersion at 150 pc scales using a new set of dust-basedα_COmeasurements and propose a newα_COprescription that accounts for CO emissivity variations across galaxies. Based on this prescription, we estimate the SFE in a sample of 65 galaxies from the PHANGS–Atacama Large Millimeter/submillimeter Array survey. We find increasing SFE toward high-surface-density regions like galaxy centers, while using a constant or metallicity-basedα_COresults in a more homogeneous SFE throughout the centers and disks. Our prescription further reveals a mean molecular gas depletion time of 700 Myr in the centers of barred galaxies, which is overall three to four times shorter than in nonbarred galaxy centers or the disks. Across the galaxy disks, the depletion time is consistently around 2–3 Gyr, regardless of the choice ofα_COprescription. All together, our results suggest that the high level of star formation activity in barred centers is not simply due to an increased amount of molecular gas, but also to an enhanced SFE compared to nonbarred centers or disk regions.